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1. Introduction 

Recent advances in the fundamentals of quantum mechanics and in quantum information 
theory pQ resulted in the renewed interest in the properties and structure of the space of 
NxN dimensional density operators (matrices). Moreover, the composite (multipartite) 



5h ' systems exhibit the effect of entanglement which still makes things more interesting 



In order to facilitate further discussion, we briefly summarize the essential properties 
of iV x N dimensional density matrices. The literature on the subject is so large, that we 
quote only two books JHE]. In the typical quantum-mechanical context the N— level 
system is associated with the Hilbert space TCn which is iV— dimensional. Quantum- 
mechanical observables are self-adjoint operators in Banach space £>(7Yat) and may be 
represented by space M.nxn of N x N hermitian matrices. Density operators p for 
iV— level systems are a subset of B(H.n) which may be denoted as B p (Hn)- A density 
operator p 6 B p {Hj^) possesses three fundamental properties 

(i) p = p\ - hermiticity; (la) 

(ii) Tr{p} = 1, -normalization; (16) 
(hi) Xj > 0, - positivity. (lc) 

where Xj - eigenvalues of density operator p. Strictly speaking we should say that density 
operators must be positive semidefinite. However, the phrase "positive" is shorter, so 
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we will use it, keeping in mind the strict sense. The set B p (TCj\f) is convex: that is, if 
Pi, Pat B P (H N ), then 

p = spi + (l-s)p 2 eB p (H N ), for se(0, 1). (2) 
Moreover, density matrices satisfy the inequalities 

Tr{p fc }<l, Tr{p fc }>^ZT, (3) 

for k > 2, while for k — 1 relation (jTTij) holds. For pure states p = the 
first inequality becomes equality (the second one is obviously satisfied). Pure states 
are also extremal - they can not be given as a nontrivial convex combination of two 
other matrices. For maximally mixed states one can find a representation in which the 
corresponding matrix is diagonal with all nonzero elements equal to 1/N. In this case 
one has Tr{p fc } = 1/N k ~ l . Finally it may be worth noting that for two-level system 
condition (fTcj) and Tr{p 2 } < 1, are equivalent. But it is not the case for systems with 
higher dimensions, that is for iV > 3. This is due to additional conditions which follow 
from requirement (JTcj) . 

Although the outlined fundamental properties of density matrices are simple and 
well-known, not much is known about the structure of the set B p (Hn). 

Only the case of iV = 2 (also called a qubit) seems to be an exception. One can 
find a one-to-one correspondence between all two dimensional p's and the set od 3- 
dimensional real vectors which are called Bloch vectors. To find such a correspondence 
one usually uses standard Pauli matrices er = (ui, cr 2 , cr 3 ) which, together with a 2 x 2 
unit matrix 1 2 , form an orthogonal basis in B{7i^) = M. 2X 2- Then one takes 



5<i> + 6-*) = 5 



/ \ 



I + 63 bi — ib 
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1 bi + ib 2 1 - 63 / 



(4) 



where b = (bi, b 2 , 63) G M 3 , is a Bloch vector and b • a = bi&x + b 2 o 2 + Density 
matrix (0J) is clearly hermitian and properly normalized. Requirement of positivity is 
equivalent to the condition Tr{p 2 } < 1 which yields the inequality which must be 
satisfied by length of the Bloch vector 



bl + bl + bl < 1. (5) 

Thus the set of 2 x 2 density matrices coincides with the Bloch ball {b e R 3 : |b| < 1} 
of unit radius. Equality occurs only for pure states which, therefore lie on the Bloch 
sphere {b G I 3 : |b| = 1}. The sphere is a boundary of the ball and the pure states are 
extremal. For maximally mixed states one has | b | = 0. 

On the other hand, even for iV = 3 the situation is not that simple. In two recent 
papers by Kimura and by Byrd and Khaneja [H] the SU(N) generators are used to 
expand an arbitrary N x N density operator. Such generators correspond to the Pauli 
and Gell-Mann matrices for N = 2, 3, respectively. The methods of construction of such 
matrices is outlined in [2] (see also |H1E3])- The expansion coefficients form a generalized 
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Bloch vector which consists of iV 2 — 1 real parameters (clearly, such a procedure is a 
generalization of the N = 2 case, because Pauli matrices are SU(2) generators). Then a 
one-to-one correspondence between density operators and the allowed generalized Bloch 
vectors is established. Moreover, it is shown that due to the requirement Tr{p 2 } < 1, the 
generalized Bloch vectors lie within a certain hyperball with a finite radius. However, 
it is also known pQ that such a hyperball contains Bloch vectors corresponding to non- 
positive matrices. This clearly indicates that the structure of the set of all allowed 
Bloch vectors (and therefore of all N x N density matrices) is not simple, not to say 
quite complicated. Investigations of the geometry of the space of density matrices are, 
therefore, difficult and complex. An example of such studies can, for instance, be found 
in @]. 

The main aim of this work is to propose a new parametrization of the set of N x N 
dimensional density matrices. Our work is somewhat similar in its spirit to the papers 
although we use other operators as a basis in B(Hn)- We will try to argue that 
the proposed representation might be useful in future applications 

The following section is devoted to the brief outline of the results obtained via the 
standard SU(N) generators by Kimura [2] and by Byrd with Khaneja |3j. We review 
these results in order to compare them to the ones presented in this work, and to argue 
that that the latter ones may be advantageous. 

In the third section we recall the concept of polarization operators and summarize 
their properties. We follow the terminology and notation used in the handbook [7j. 
Similar operators (although in a less transparent notation) are also discussed in the 
classic book by Biederharn and Louck jH]. These authors also outline the application 
of the polarization operators to the expansion of density operators. We elaborate on 
the idea in the fourth section. Moreover, we discuss some properties of the expansion 
which seem not to be documented in the literature. Namely, we study the analog of the 
generalized Bloch vector in the light of the conditions imposed upon a density operator. 

Investigations of positivity of the density operator are much more involved than 
checking hermiticity or normalization. Hence, the fourth section is devoted to this 
issue. We derive general expressions which allow construction of positivity conditions 
for any N. 

In the two next sections we employ the developed procedure for a qubit (N = 2) and 
for a qutrit (N = 3). The first case is simple, while in the second one we show that the 
positivity requirements are quite restrictive. We study two-dimensional cross sections 
of the space of generalized Bloch vectors which possess a complicated and asymmetric 
structure. 

Finally, in the last section we give some concluding remarks and indicate some 
possible future applications and further developments. 
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2. Standard SU(N) representation 

The authors of the papers [2J dj present very similar ideas using, however, somewhat 
different notation. Their results may be viewed as a generalization of the N = 2 case 
to higher dimensions. The idea is to expand the density operator in some, suitably 
chosen basis of orthogonal (and traceless) matrices. Such a basis in B(Hn) is given 
by the set of matrices which are the standard generators of the SU(N) group together 
with the unit matrix ljy. The SU(N) generators have the following properties (for 
i,Jfe = l,2, ...,iV 2 -l) 

(i) Ak = A k , - hermiticity; (6a) 

(ii) Tr | A k j = 0, - tracelessness; (66) 

(iii) Tr | Aj Afc| = 2 <5jfc, -orthogonality. (6c) 

Property (j6j| entails easy normalization of the density operator, while (iii) (being the 
Hilbert-Schmidt scalar product) ensures that the given matrices indeed form orthogonal 
basis in B(TLn)- 

The commutators and anticommutators are given as (summation rule holds): 

[Aj, A k ] =2if jkm A m , (7a) 
4(5-, A 

[A;, A fc ] + = — p 1 N + 2g jkm A m , (76) 

where fjkm is a completely antisymmetric tensor, while gjkm is a completely symmetric 
tensor. As mentioned, construction methods [2JIH1EII| of matrices Aj-SU(N) generators 
and of the structure constants fjkm, gjkm are known, but fairly complicated. 

Then, following Kimura |2. one uses operators Aj in B(Hn)) one can write for 
the density operator (Byrd and Khaneja [3 adopt slightly different normalization of bj 
coefficients): 

N 2 -l 

where bj G R, give the generalized Bloch vector b = (pi, 6 2 , , &iv 2 -i) _1 - Due 

to the properties ()6aj) and (|6*c|) one sees that this expansion clearly yields a hermitian 
and normalized matrix. Moreover, requirement Tr { p 2 } < 1 implies that the vector b 
must satisfy the condition 



\ 



N 2 -l 



So we see that vector b characterizing an arbitrary A^-dimensional density operator 
must lie within (A^ 2 — l)-dimensional sphere specified by the inequality Q. 

As already mentioned above, for cases where N > 3 requirement © isn't equivalent 
to Aj > 0. Positivity imposes some additional restrictions on vector b. Due to that, 
generalized Bloch vectors constitute a subset inside the hypersphere. This set has 
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complicated and asymmetric structure as briefly discussed in 0. Positivity of density 
operator is then investigated along the same lines as we are employing in this work. 
Hence, we postpone the discussion of this issue to subsequent sections. 

3. Alternative representation for Bloch vector 

3.1. Polarization operators 

In this section we recall the facts given in the handbook |2j (chapter 2.4.) on quantum 
theory of angular momentum. For sake of completeness of this paper we define the 
concepts, introduce notation and briefly give some additional comments avoiding further 
reference to the given source. 

Let £j denote, for a given but fixed value of j, the space spanned by eigenvectors 
\jm) of angular momentum operator. Number j is half-integer or integer, so it can 
take values: j = ~, 1, §, 2, . . ., while obviously m = j, j — 1, . . . , —j + 1, — j. Since for a 
given j there are 2j + 1 vectors \jm), the space Sj is (2j + l)-dimensional. Vector \jm) 
is then represented by a column vector of 2j + 1 components, 2j of them being zeroes 
while the m-th one (m = j, j — 1, . . . , — j, in that order) is equal to 1. From now on 
we will identify £ 3 ■ = Hn = T~C2j+i, thus specifying Hilbert space for a N = 2j + 1-level 
system. 

In the operator space B(H2j+i) we introduce the following set of operators: 

Tlm(j) = M±± J2 C ^,LM\jm)(jm'\, (10) 

V m,m' 

which are called polarization operators. Cj^' lm are Clebsch-Gordan coefficients 
(CGC). Due to the properties of CGC one immediately sees that numbers L and M are 
always integers and take the following values: 

£ = 0,1,2, ,2j, M = -L, L, (11) 

It is straightforward to see that there are 

Y%Lo( 2L + 1 ) = ( 2 J + !) 2 polarization 
operators. This suggests that these operators constitute the basis in B(Hn)- After 
discussing the properties of polarization operators, we will show that this is indeed the 
case. We also note that CGC are well-known, fully documented and easily computed 
by the computer. We proceed with listing the fundamental properties of the introduced 
polarization operators. 

Employing the symmetry properties of CGC one easily shows that the hermitian 
conjugate of operator Tlm(J) ^ s gi yen as 

TluU) = (-!) M Tl-uU), (12) 
so, in general, polarization operators are nonhermitian. As we will see later this does 
not pose any serious difficulties. On the other hand, operators T L0 (j) are diagonal and 
hermitian. 
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For any j > | and L = 0, M = 0, one has 

T oo(j) = ^ 2 j + 1 ^ JV=2 ^+ 1 ' ( 13 ) 

so operator T o(j) is always proportional to identity while multiplication coefficient 
depends on j. 

Operators Tlm are traceless, in the sense that 

Tr{ Tlm U)} = v^j + l hohio- (14) 
Indeed, directly from the definition (JTU|) we have 



MWi)} = E C Zlm = f^^q- L ,(15) 

which follows from the condition M + m = m which must be satisfied by nonvanishing 
CGC. Performing the last sum one obtains relation (|14|). 

Computation of the trace of two polarization operators is more tedious. Using the 
definition (fTUj) and employing the known expressions for the sums of products of CGC, 
one gets 

Tr{T LlMl (j)^L 2 M 2 (j)} = (-l) Ml 6 LlL2 S Mi -m 2 (16) 

From this relation and due the property ()12|1 the Hilbert-Schmidt product of polarization 
operators is given as 

T L 2 M 2 (j)} = S LlL2 5 Ml M 2 ( 17 ) 

From this we conclude (as it was also in the case of standard SU(N) representation) 
that operators Tlm form an orthonormal basis in space B(HN=2j+i)- 

Similar calculations can be performed to obtain the traces of multiple products of 
polarization operators. Then for any n one finds 

Tr{ T LlMl (j) T L2 m 2 (j) T LnA / n (j) } = 

— iil2...n ^ LiM-i, jrn ^ L 2 M 2 , jm+m ■ ' • °i n _iJW„_i, jm+At„_ 2 ^ L„M n , jm+Hn-!' K^ou) 

m 

where we have introduced the following notation 



^ = f Ll + 1 \ 2 r + ^ K + 1)) and W = 5>- (186) 

It is perhaps worth noting that the properties of CGC imply that Ym=i ^ = ^' ^ 
otherwise then the trace is zero. It can be seen from the last of CGCs in (118 oj) where 
M n + \i n -\ — 0- Certainly, relations (fHj) and (fTEJ) follow from the general expressions 
(|3.1)1 . The traces of products of polarization operators are expressed by a relatively 
simple formula. This should be confronted with results of Byrd and Khaneja P] who 
find quite involved expressions for the traces of the products of A$ matrices up to n = 9. 
In the present case, formula ()18aj) is valid for any n. 
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The product of two polarization operators is computed in the similar manner. 
Adapting the known sums of CGC to the present needs, we express the product by 
a combination of polarization operators 



(J)TWJ) = 5>l) W V(2£i + 1)(2L 2 + 1) ' Ll L " U 
La 13 3 3 



x C LimIl 2 m 2 TlsMsO')) (19) 
which involves Racah 6j-coefficient. This relation can also be used to find trace of the 
product of two polarization operators. Since TL 3 M 3 (j) are traceless for L 3 > 0, only 
L 3 = M 3 = contribute to the sum. Then C°° M l 2 m 2 i m P nes that M x = — M 2 and 
the properties of the remaining coefficients yield relation (fTo^l . Equartion (jl9j) allows 
extension to multiplication of any number of TlmU) operators in an easy but somewhat 
arduous way. 

Another essential properties of any set of operators are given by their commutators 
and anticommutators. They follow immediately from expression (fTT^) and from 
symmetry properties of CGC, and are given as: 

[TWJ), UOlr v / (2^i + l)(2L 2 + l) J2 (-1) W3 




x [ 1 ± (-1)^+* ){ - \ Cgjg, L2M2 T /..',/,(./). (20) 



These relations also show that polarization operators (slightly differently normalized) are 
just another realization of the representation of the SU(N) group. The above relations 
specify structural constants other that those used in j2HH]- This explains why we have 
called the representation used in these papers "a standard one". It is worth stressing 
that properties of TlmU) operators follow directly from the symmetry properties of 
Clebsch-Gordan coefficients and from properties of other quantities well-known from 
the impressive literature on quantum theory of angular momentum. We also note that 
Biederharn and Louck give some interesting comments on the structure of the present 
alternative representation of SU(N) group generated by operators TlmU)- Discussion 
of these issues, although quite interesting, goes beyond the scope of this work, but 
perhaps deserves some further study. 

Standard generators of SU(N) formed a basis in B(TIn). The same applies 
to polarization operators, therefore they can be clearly used to expand operators. 
Let A G B(TCN=2j+i) with number j specifying the dimension of the Hilbert space: 
N = 2j + 1. Since j is fixed for a given physical system, henceforward we will not write 
it explicitly where it is not essential. 

Obviously, due to completeness of states \jm) which span the space Ti,N=2j+i we 
can decompose any operator A acting on this space as follows 

A = ^2^2A mm >\jm)(jm'\, (21) 
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with A mm i = ( jm | A | jm' ) being the corresponding matrix element of A. The set of 
matrix elements A mm i completely describes operator A in a given basis. 

Since polarization operators Tlm constitute another orthonormal basis in B(TCn) 
we can decompose any operator in this basis, thus writing 

2j L 

a = Yl Alm Tlm ' ( 22 ) 

L=0 M=-L 

Due to orthonormality relation (|17|) the components Alm of the above expansion are 
given as 

A lm = Tt{V lm A}. (23) 

These coefficients completely describe operator A in the basis of Tlm operators. 

Employing the well-known orthogonality relations for CGC one can find the 
relations between A mm i and Alm 




2L + 1 

2j _|_ I ^jm', LM Amm>, (24a) 



L=0 M=-L 



or _|_ i _ 

C]Z,lmAlm. (246) 



3.2. Decomposition of density matrix 

Clearly, a decomposition such as discussed above may also be applied to density operator 
describing the N = 2j + 1-level system (see also chapter 7 of Reference jHj). Since 
operator T 00 plays special role (due to the relation ([13])) we will write it outside the 
sum. Hence we have, for the considered density operator, the following expansion 

2j M=L 

P = V 00 Too + £ Yl Vlm 1 lm = V °° T °° + ^ ' T " ( 25 ) 

L=l M=-L 

where V is again called a generalized Bloch vector, similarly as it was done in case 
of the expansion in terms of standard SU(N) generators. This vector consists of 
iV 2 — 1 = (2j + l) 2 — 1 components ordered as 

V = (Vi-i, V 10 , V n , V 2 - 2 , V 2 -x..., Vlm). (26) 

The component Vlm of the generalized Bloch vector is given similarly as in (}2*3*j) . that is 

Vlm — Tr{ T LM p}. (27) 

Polarization operators are traceless. This fact, together with the normalization condition 
dljj and with relation (JHJ), imply that 

1 = Tr{ p } = V 00 Tr{ T 00 } = V 00 y/2j + l, (28) 

and therefore, for any fixed N = 2j + 1 we have 

Voo = ; 1 = -!=, (29) 
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which is then also fixed. Hence the decomposition (|25|) may be rewritten as 

p = ±i N=2j+1+ V.f } (30) 

where the product V • T is specified by (ffijj) . Polarization operators are traceless, thus, 
the required normalization of the density operator is always preserved. Moreover, this 
explains why we have restricted the generalized Bloch vector V to L > 1, as indicated 
in definition ()26|) . 

Density operator must be hermitian, while polarization operators Tlm are not. 
Therefore, components Vlm of the generalized Bloch vector V are in general complex. 
Due to relation (j!2j) the hermiticity of the density operator implies that the complex 
conjugates of components Vlm are given as 

V*LM = (-l) M V L - M (31) 

Components of the generalized Bloch vector have the following properties: 

• there are N — 1 = 2j components Vlo (with L > 1). Since T^ = T lo numbers Vlo 
are real; 

• there are N(N — 1) components Vlm (with L > 1 and M ^ 0). Relation (j31j) 
implies that: (i) for even M, if Vlm = ct + i (3 (with a,/?6l) then Vl-m = at — i/3; 
(ii) for odd M, if Vlm = a + i j3 (with a, /3 6 K) then Vl_a/ = —a + i /3. So these 
components are fully represented by N(N — 1) real numbers. 

Therefore, we conclude that relations (|31|) ensure hermiticity of the density operator and 
that V is given by iV 2 — 1 real numbers, just like the corresponding generalized Bloch 
vector in standard SU(N) representation. Thus, we can say that using new type of 
operator basis we obtain alternative decomposition of iV-dimensional density operator, 
which preserves two requirements: hermiticity and normalization. The question of 
positivity is much more involved, so we shall discuss it separately in the next section. 

Before doing so, we would like to add some additional remarks. The fact that Vlm 
are complex has another interesting consequence, namely, their values as defined in (|2"Tj) 
cannot be measured directly. We have to construct a set of hermitian observables using 
polarization operators. The solution is simple, for instance, we can take the following 
hermitian combinations: 

Qlm = Tlm + T^ M , (32) 
Qlm = ^(Tlm — T iM ), 

and now Vlm are expressed by the expectation values of these observables 

Vlm = 1(Qlm) + ^(Qlm), (33) 

which, in particular, yields Vlo — (Tlo)- Specific experimental needs may require 
construction of yet anotherset of hermitian observables, however, we have shown that 
this is clearly possible to do so. Hence we see that the fact that Tlm are nonhermitian 
does not present any difficulty 
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Moreover, for the sake of future needs, we define the scalar product of two 
generalized Bloch vectors Vi and V2 by the following formula 

Vx-V 2 = £ E V&(vffi). = {-1) M VHV?} M (34) 

L=l M=—L L=l M=-L 

As a next remark we give a restriction on the length |V| 2 = V • V (in the sense 
of (J33j) )- Any density operator must satisfy inequalities (J3J). Employing expansion (JHOJ) , 
using the tracelessness of polarization operators and expression (|TB|) we arrive at the 
equation 



Tr{p 2 } = l + Tr|(v-f) 2 } 



2j L 



A? + E E WZm = ^ + |V| 2 < 1. (35) 

L=l M=-L 

Note that we have Tr <| ( V ■ ]> = | V| 2 . Relation (J3l)j) yields the mentioned restriction 

(36) 




2j + V 

This result is analogous to relation (jHJ) obtained in the case of the standard SU(N) 
representation. Hence, we can say that that all generalized Bloch vectors V form a 
subset within a hypersphere of radius R = y(N — 1)/N. The structure of this subset 
(for N >3) is, as we know, pretty complicated because requirement of positivity imposes 
additional restrictions on this subset. Nevertheless, pure states (for which Tr{ p 2 } = 1) 
lie on the surface of the hypersphere of radius R. On the other hand, maximally mixed 
states (for which Tr{ p 2 } = l/N) correspond to |V| = 0. We can say that the shorter is 
vector V, the density operator represented by it corresponds to a "more mixed" state. 
The length of V can be, thus, considered as a kind of measure of " mixedness" jH] . 

Finally, we note that for two density operators p\ and p2 one has < Trjpi^} < 1 
(see 0j). Using two expansions (jSOjl and due to tracelessness of T £M 's we obtain 
inequalities for corresponding vectors Vi and V2, namely 
1 - N -I 

-7r* v, - v> *— ' (37) 

If both these vectors describe pure states then their lengths are equal ■sj (N — 1)/N and 
the above inequalities are expressed in terms of cos a - cosine of the angle between two 
considered vectors. Thus, we obtain for this case 

- — - — < cos a < 1. (38) 
N- 1 - - v ' 

This result reproduces the ones obtained in [2] and 0. It confirms the notion that the 
set consisting of all allowed generalized Bloch vectors is complicated and quite likely 
asymmetric. 
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4. Positivity of density operator 

Any density operator, apart from being normalized and hermitian, must also be positive. 
Requirement of positivity may be formulated in many ways. For example, operator p is 
positive if 

V|^)gHat => (ip\p\tp) >0. (39) 

Another formulation is given via the eigenvalues, as it was stated in (fTcj). We shall 
concentrate on the latter approach since it was also used by the authors of papers [21 E]- 
In this manner we would be able to compare our results with those of Kimura and Byrd 
with Khaneja. 

Thus, we need necessary tools to investigate the eigenvalues of the iV-dimensional 
density operator. Such a tool is provided by the characteristic polynomial of a variable 
A G 1: W(X) = det(Aljv — f>) (see Ref.jB]). Note that, in comparison with the usual 
notation, we have changed the sign. It can be shown that the polynomial W(X) may be 
written as 

N 

W(X) = Y,(-V k Sl N) X N - k , (40) 

with Sq N ^ = 1. It is also worth noting that = det p. Coefficients S[ , (k = 

1, 2, . . . , N), are constructed recursively by the Newton's formula 

kSi N) = ^(-ir- 1 S^ m Tr{p m }. (41) 



m=l 



Obviously S{ = Tr{p} = 1 due to normalization of the density operator. Computation 
of subsequent s is straightforward. Several initial quantities S% are as follows 

= \~\ Tr{p 2 }, (42a) 

S f ] = \-\ Tr{p 2 } + ^ Tr{p 3 }, (426) 

S[ N) = Y A ~\ Tr{p 2 } + I (Tr{p 2 }) 2 + \ Tr{p 3 } - \ Tr{p 4 }, (42c) 

and will be useful in the study of some examples in subsequent sections. The same results 
(although in different notation) are also given in Refs.j3 El- Having constructed the 
coefficients Sk of the characteristic polynomial (}4*0|) of the density operator p 6 B p (T-Cn), 
we can address the question of positivity. The answer is supplied by the following 
theorem 

{S ( k N) > 0, for all k = 1, 2, . . . , N } ^ {V,- Xj > 0, that is p - positive}. (43) 

Hence, to check whether a given operator p is indeed positive, one needs to check the 
positivity of the corresponding coefficients . On the other hand, requirement that 
, (k = 2,...,N) are nonnegative imposes restrictions on the components of the 
vector V, thereby inducing a complex structure on the set of all allowed V's. 
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Certainly the condition > (see ()42ajl ) implies Tr{p 2 } < 1 and therefore 

reproduces the requirement < |V| 2 < (JV — 1)/N, as already discussed in (JHHJ). So, 
for a qubit (when j = 1/2 so that N = 2) the latter requirement is indeed equivalent 
to the requirement of positivity. For higher dimensions (k = 2, . . . , N), and then 
the condition imposed upon Tr{p 2 } is necessary but not sufficient to ensure positivity. 
The higher (that is for k > 3) must also be checked for nonnegativity. This simple 
remark, strangely enough, seems not to be noticed in papers and [3]. 

Quantities are easily computed provided the traces Tr{p fc } are known. Density 
operator is normalized and equation ()35|) gives Tr{p 2 } = 1/N + |V| 2 . For k > 3 we 
employ expansion ()30|) which yields 

i :•■ i i I ; -i v ■ v-i) ; no 




Using Newton's binomial and tracelessness of polarization operators we have 
,k^ f k \ T m 



Tr{p fc } 



-m 



m=0 



_J_ fc(fc-l) |y|2 + W k \ _Tm_ (45) 



m=3 



where we have denoted 

T m = Tr{(V-T) L (hi) 



which, due to hermiticity of p are real. In expression (J45j) we understand that 
(V • T)° = 1 7v which entails that T = N. Moreover, tracelessness of polarization 
operators implies that T\ = and relation (|35|) gives T<i = |V| 2 . So we can say that the 
problem is now reduced to computation of the quantities T m for m > 3. Directly from 
the definition, it follows that one can write 

2j Li 2j L k 

T *=E E - E E ^ 1 M 1 ---Vl fcA4 Tr{T LlMl ...T LfcM J, (47) 

L 1= lAfi=-Li L k =l M k =-L k 

and since the multiple trace is known (see formulas ()3.1)) ) we can find any and 
therefore the traces Tr{p fc }. Resulting expressions are complicated but the multiple 
trace is not zero only when X^i=i % = which greatly reduces the number of terms. 

Before constructing explicit expressions for coefficients we write down traces 
of p k for k — 3, 4. They will be useful later and are as follows 

Tr{p 3 } = ^ + !|V| 2 + T 3 , (48a) 

Tr{p 4 }=^ + ^|V| 2 + |T 3 + T, (485) 

Coefficients follow by combining the recurrence relation (j4*Tj) with the first of 

equations (j43|) 

^^EHi'-sSEi.)^ («) 

fc=l m=0 
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Writing out the k = 1 term, we obtain 



n 



c{N) _ n{N) 

J n — °n-l / j 



m=0 



1 

m 



T 

' it 



N 1 - 



n K 



k 
m 



T 

J. 77 



k=2 



m=0 



(50) 



Certainly, the second term contributes only for n > 2. Since Sq ^ — S± 1 — 1 we can 



,(TV) 



safely assume that n > 2. Then, we note that T 



n 



n 

sw = sS9 1 + Y(-i) i - 1 s 



k=2 



(TV) 
n—k 



N and T x 

k 



0, hence 



k 
m 



(51) 



which is the sought recurrence relation for coefficients S^K The traces T m can be 



computed as discussed above. The first nontrivial coefficients are (with T2 = |V| 2 )) 

T,, (52 a) 



S. 



(TV) 



iV-1 



s. 



(TV) 



2N 2 
(N — 1)(N — 2) 



S. 



(TV) 



6N 2 

(N - 1){N - 2)(N - 3) 



N-2 1 
"2AT T2+ 3 T3 ' 



(N-2)(N-3) 



2AN 3 



AN 2 



To 



N-3 m 
3A^ 3 



(526) 



(52 c) 



As the theorem states, positivity of p is equivalent to the conditions that > 
for all k — 1, 2, . . . , N. Relation (|51J1 allows one to compute these quantities for any finite 
N. These computations might be lengthy or tedious, but otherwise straightforward. 
This follows from expression (|%7j) which together with relations (|3.1|) allow us to compute 
traces T^. Thus, we conclude, that the proposed approach to the parametrization of 
N- dimensional density operator can be applied in a closed form for any N. On the other 
hand, Kimura [2] gives specific expressions only for N < A while Byrd and Khaneja |3j 
give up to iV < 9. Our presentation is free from such restrictions. We give quite specific 
expressions valid for any N. 



5. Example: Qubit 

The formalism introduced above can now be applied in some specific cases. The simplest 
one is a qubit and it corresponds to j = 1/2, hence to iV = 2j + 1 = 2, which within 
the "standard SU(N) framework was described by Pauli matrices and Bloch vector. In 
this case prescription (fTUj) yields L = 0, 1 and M = —L, . . . , L for each L. The matrices 
representing polarization operators are of the form 



Is 1 / 1 \ _ ,1 





fl 




V° 1 




fl 




\ - 



T,+i( 
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Then, according to (|3U|) we can represent the 2- dimensional density matrix by 3-dimen 
sional vector V, 

\ 



P 



1M 1 1M 



( 1^1 

2 + V2 



M=-l 



(54) 



Denoting Vi 
parameters 



-V u 

1 _Vn 
V 2 ^ / 

ieK, Vii = a + if3 G C and using relation we get in terms of real 



— V* 



/ 1 

2 

V " 



x 

"7! 

a + i(3 



— a — i/3 



\ 



1 x 

2 + 7! / 



(55) 



(2) 

this matrix is clearly hermitian and normalized. The requirement of positivity £2 > 
is equivalent to the condition Tr{p 2 } < 1, which then gives 

1 



V 



Vf + 2\V n \ =x 2 + 2(a 2 + (3 2 ) < 



2' 



(56) 



which is an analog of relation (jHJ). The surface x 2 + 2(a 2 + f3 2 ) = 1/2 can be represented 
parametrically 



Re{^n} 



sint 



cosu, Im{Vii} 



sint 



sm u, 



cost 



10 



(57) 



where t G [0, 7r] and u G [0, 27r]. This is a prolate spheroid which in the present case 
corresponds to the "standard" Bloch sphere. All allowed vectors V representing a 2- 
dimensional density matrix lie within this spheroid, while pure states occupy its surface. 
Hence, the proposed description of the density matrix yields results fully equivalent to 
the "standard" one. 



6. Example: Qutrit 



6.1. Construction 



In this section we discuss the next example. We investigate the 3-level system, sometimes 
called a qutrit. For such a system we have to N = 3, and that entails j = 1 in the spirit 
of section 3. The operator basis is spanned by 9 polarization operators Tlm{1) with 
L — 0, 1, 2 and M = —L, . . . , L Applying the rule (fT0|) we construct the corresponding 
matrices. They are as follows 



TV 



1 



1 oo : 

(° 



V° 



i 

7! 



(i 

o 







1 



10 



1 

7! 



/ 1 

o 





- 1 / 



1-1 



1 

7! 



/o 
1 



\ 



1 0y 
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1 2+2 



/ 1 \ 






T2+1 — — ^= 



1 

3 



(0 



v° 



-1 o\ 

1 
j 



T20 — 



1 






\ 
-2 
l) 



I 2-1 



1 

72 



/ 















/ 











\ 




1 










, T2-2 = 















V 





-1 





J 




\ 


1 








/ 



(58) 



Thus, as expected, matrix T o is proportional to the identity one, and the corresponding 
density operator is written, according to the prescriptions (I2HJ) and (jSHJ), as 

1 2 L 1 

P = ^^ + Y1 J2 VlmT lm = -1 3 + V-T, (59) 



where V = (Vi_i, V w , Vu, V2-2, V2-1, V20, V21, V22) (ordered as written) is a generalized 
Bloch vector representing a qutrit. Moreover, we note that according to relation (f2*H|) 
we can also write V 00 = l/v3. Eight complex components Vlm must satisfy relations 
(|3ip. As it follows from the comments given after these relations they are also specified 
by 8 real numbers. Since dealing with real quantities is simpler and more transparent, 
we introduce the following notation 



V10 
V20 



x, 

y, 



v n 
v 2 i 



—V* 

v i-i 

—Vo 



a + ib, 

, 2 _ x _ ai + ifc, V 22 = V2-2 = «2 + i#2, (60) 

with x,a,b,y,ai, Pi,a2, P2 G With the aid of this notation we can explicitly write 
down the the density matrix for a qutrit 

/ 



I x y 
3 v/2 v 7 ^ 

2 — ib oti — ifli 



a + ib ai + if3\ 



\ 



\ 



y/2 V2 
a 2 - ifa 



V2 

1 

3 



y/2 



a 2 + i{3 2 



a + ib ct\ + if3\ 



V2 



V2 



ib ct\ — i/3i 



1 x y 
3 ^2 



(61) 



/ 



y/2 V2 

which is clearly hermitian and normalized, as it should be. 

To identify matrix (JBTj) as a true density matrix we must be sure that it is positive. 
The positivity conditions correspond to the inequalities > and > 0. The first 
one, as already discussed, is equivalent to 

,2 



< IVI 



2 

< 3' 



(62) 



with the left inequality being trivial due to the definition of |V| 2 . 
immediately from (J52 6)1 and it reads 

3 9 2 1 1 



follows 



S. 



(3) 



+ T 3 



> 0. 



(63) 
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The next quantity necessary to investigate the positivity of the density matrix (J6*T|) is 
T 3 = Tr{(V • T) 3 }. It follows from relations (|47|) and (|18a|) . Computation of the triple 
trace is a bit tedious but otherwise straightforward since the CGC are tabulated in [Jj. 
So, with identifications (JBTJj) we obtain 



one expresses S3 via the introduced 8 real variables. Moreover, one easily checks that 
= det p, as is should be. 

6.2. Parametrization with two nonzero variables 

General analytical discussion of positivity conditions and together with (|54*|) 
and ()65|) seems to be extremely difficult if not virtually impossible, because there are 8 
real parameters. Therefore, we will restrict our attention to a simpler case. Namely, we 
will assume that only two of real parameters (x, y, a, b, ax, Pi, 0.2, P2) are nonzero while 
the other six ones are put to zero. Similar procedure was employed by Kimura [2j and 
Byrd with Khaneja [H|- Then, there are 28 different pairs of nonzero parameters. We 
shall show that these 28 pairs split into 7 distinct types. 

We shall denote the pair of nonzero parameters by (s, t) and next we will indicate 
to which pairs taken from the set (x, y, a, b, ax, Pxi a 2, P2) it corresponds. Then we will 
give then conditions (}6*2*j) and (|63j) written in terms of parameters s, t. The formal 
expressions for these conditions are common to all representatives of the given type. To 
discuss (J53j) we introduce the quantity 



which will be a function of two variables s and t and the inequality obviously follows 
from (|63jl . After presenting the basic properties of each of the seven types we will briefly 
state the properties which are common to all of them. 

6.2.1. Type I Type I corresponds to only one pair of parameters, that is to (s,t) = 
(x,y). Then, requirements (J62j) and (J66j) translate into the following ones 



3x 2 y — y 3 

7$ 




a 2 _ h 2 _ a 2 _ p2) 




1 1 , o 

F — V 2 + T 3 > 

9 2 1 1 3 ~ 



(66) 




x 2 + y 2 < -, 



(67a) 



F(x,y) 



1 x 2 + y 2 3x 2 y — y 3 



> 0. 



(675) 



9 2 v/6 




Figure 1. Type I (a) left and type II (b ) right. Allowed states lie on the solid lines 
and inside the triangles. Small circles - pure states. The grid unit corresponds to 
actual value equal to 0.1. 



Requirement (j67oj) restricts allowed values of x, y to the circle of the radius R = a/2/3 
(dashed line). Then (|67fr|) implies that the allowed points must lie within or on the 
triangle indicated by solid lines in Figure 1. We also note that for this type there are 
three possible pure states (±1/a/2, 1/VQ) and (0, —a/2/3). They are indicated by small 
circles in Figure 1. 

6.2.2. Type II Type II has two representatives (s, t) = (y,a 2 ), (y,^)- Then relations 
(1521 and © give 

2 

|V| 2 = s 2 + It 2 < -, (68a) 
3 

^>4-T- i2 + ^ £a (686) 

Inequality ()68aj) places the allowed values of parameters on and inside the ellipse with 



semi-axes of lengths a/2/3 and a/1/3 (dashed line in Figure 1 (right)). Inequality ()68&j) 
restricts s and t to the triangle drawn with solid lines in Figure 1 (right). There are 
also three pure states (\/l/6, ±a/1/2) and (—a/2/3, 0). 

6.2.3. Type III Type III is specified by two cases (s,t) = (a, a 2 ), (/?i,«2)- Then from 
(JH2J) and flESJ we get 

\V\ 2 = 2{s 2 + t 2 )<\, (69a) 

F(s,t) = - - s 2 -t 2 + 3s 2 t > 0. (695) 
9 

The first condition gives allowed values of parameters s and t within a circle of radius 
R = a/1/3 (dashed line in Figure 2 (left). The requirement (|69fo|) restricts the values of 
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-7 -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 -7 -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 



Figure 2. Type III (a) left and type IV (b) right. Allowed states lie on the solid lines 
and inside the region restricted by a horizontal straight lines and parabolas. Small 
circles - pure states. The grid unit, as previously, corresponds to actual value equal to 
0.1 



s, t to a region below the straight line t = 1/3 and above the parabola t = 3s 2 — 1/3. 
This region is indicated by a solid line in Figure 2 (left). This type allows for two pure 
states (±v2/3, 1/3) which are denoted by small circles. 

6.2.4- Type IV Type IV is similar to the previous one. In this case we also have two 
cases (s,t) = (6,02), (011,012). Then from (|62|) and (|66|) it follows that 

|V| 2 = 2(s 2 + t 2 ) < -, (70a) 
3 

F(s,t) = - - s 2 -t 2 -3s 2 t > 0. (706) 
9 

The imposed conditions are thus similar. Only now, inequality (j706|) implies that 
the allowed points are above the straight line t = —1/3 and below the parabola 
t = —3s 2 + 1/3 drawn as solid lines in Figure 2 (right). Two pure states correspond to 
points (±V2/3, -1/3). 

6.2.5. Type V Type V has four representatives, namely (s, t) = (y,a),(y,b),(y,azi), 
(y,Pi). Next, from (|fi2~j) with (jHEJ) we have 

|V| 2 = s 2 + 2t 2 < -, (71a) 
o 

F(s, t) = ----t 2 - 3t2s + s3 > 0. (716) 
1 ; 9 2 VQ ~ 

The first requirement puts the allowed points within an ellipse with semi-axes of lengths 

a/2/3 and a/1/3. Condition ()716|) implies that the values of s and t are to the right 
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Figure 3. Type V. Allowed states lie on and inside the solid ellipse. Small circles - 
pure states. The grid unit equals to 0.1 



of the straight line s = —a/2/3 and within an ellipse the center of which is shifted to 
the point (— l/2y/E, 0). Semi-axes of this ellipse are equal to a/3/8 and a/1/8. So, the 
allowed points are within (and on) the smaller ellipse which is drawn with a solid line 
in Figure 3. Here, there is only one pure state at the point (—a/2/3, 0) indicated by a 
small circle. 




Figure 4. Type VI (left) and type VII (right). The allowed states lie on and within 
the solid curves. There are no pure states. The grid unit equals 0.1, as always. 



6.2.6. Type VI This type is most numerous, it is specified by 11 cases. They 
are (s,t) = (a, b), (a, ai ), (a, ft), (a, (3 2 ), (b, a ± ), (b, ft), (b, (3 2 ), K, A), (a 1} (3 2 ), (A, fo), 
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(a 2 ,/3 2 )- Then, requirements ()62jl and yield 

|V| 2 = 2(s 2 + t 2 )<^, (72a) 

F(s,t) = - - s 2 -t 2 > 0, (726) 
9 

because in this case we find T 3 = 0. Both conditions specify circles. Inequality ()72aj) 
gives a circle with radius R = a/1/3 while (J726)) yields r = a/1/3. The smaller circle 
(more restrictive) is drawn solid in Figure 4. No pure states are allowed here. 

6.2. 7. type VII The last type is given by 6 pairs, that is by (s, t) = (x, a), (x, 6), (x, at), 
(x,/3i), (x,ac2), (x,^). Then relations (J52~j) and give the requirements 

|V| 2 = s 2 + 2t 2 < \ (73a) 

F( S ,t) = i-^-t 2 >0, (736) 

since in this case we also have T 3 = 0. Condition (j73aj) restricts the allowed values of 
s and t to an ellipse with semi-axes equal to a/2/3 and a/1/3. The second condition 
restricts the allowed points to within and on an ellipse with semi-axes y/2/3 and 1/3. 
The smaller ellipse is drawn solid in Figure 4. As in the previous type there are no pure 
states. 

6.3. Features common to all seven types 

All discussed types are characterized by two-parameter behaviour of |V| 2 and F(s,t), 
the latter being proportional to det p (up to a factor of 1/3). In all considered cases 
requirement F(s, t) > is more restrictive than the condition imposed upon the length 
of the generalized Bloch vector V. That confirms the idea that the set of all density 
matrices is a proper subset of the hyperball determined by requiring that |V| 2 < 2/3. 

The allowed values of parameters lie on and inside the solid contours which 
correspond to F(s,t) = 0. Since for pure states we have |V| 2 = 2/3 and det p = 
it is not surprising that pure states are situated at the extremal points of solid contours. 

Inside these contours F(s,t) is obviously positive and for all cases attains its 
maximal value equal to 1/9 at the point (0,0) which corresponds to V = 0, that is 
to a maximally mixed state. 

7. Final remarks 

We have presented and discussed the representation of the N x N dimensional density 
matrix in terms of polarization operators TlmU)- This idea is not entirely new (see, 
for example [3 El), but the usefulness of this expansion seems to justify our recollection 
of known facts. On the other hand, we have discussed the important issue of positivity 
which seems not, in the context of polarization operators, to be considered in the 
literature known to us. 
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Usefulness of the presented approach, as it seems to us, stems mainly from the fact 
that polarization operators are expressed in terms of quantities well-known from the 
quantum theory of angular momentum. Connection with this theory greatly facilitates 
all considerations and allows derivation of expressions valid for any N = 2j + 1. For 
example, it is straightforward to find commutators, products, traces over products, etc. 
This can be done analytically, but also with the aid of computer programs allowing 
symbolic mathematics. These possibilities seem to indicate that the discussed approach 
is indeed useful in practice. 

Polarization operators TlmU) constitute a basis in the space of N x N dimensional 
operators. Hence, the density operator can be expanded as in (fHU|) . This establishes a 
relationship between N x N dimensional density matrices and V's - generalized Bloch 
vectors. Expansion (|30j) automatically ensures proper normalization of the density 
operator. This is due to the fact that polarization operators are traceless. Next, 
requirement of hermiticity implies that complex components V£ M = (— 1) m Vz,-m which 
reduces the number of independent real parameters to iV 2 — 1 and fixes the dimension 
of the space of generalized Bloch vectors. 

Density operator must be not only normalized and hermitian but also positive. 
This problem seems not to be discussed earlier a in terms of polarization operators. We 
have done that and presented the formalism allowing one to check whether the positivity 
conditions are met by a given hermitian and normalized matrix. Successive positivity 
requirements (that is inequalities S% > for k = 2,3,..., N) are derived for any 
N. These requirements specify and restrict the space of generalized Bloch vectors which 
represent true density matrices. Expressions (|47|) together with (jlSaj) allow computation 
of the quantities T m = Tr{(V- T) m } which, in turn, are used to compute coefficients 
as in (|5ip. In particular, requirement S2 > entails Tr{/3 2 } < 1 which, in terms of 
Bloch vector yields | V| 2 < (N — 1)/N. So, all vectors V lie within a N 2 — 1 dimensional 
hyperball of radius equal to y/(N — 1)/N. However, further conditions > (for 
k > 3) severely restrict the set of allowed Bloch vectors. This set is a proper subset 
of the mentioned hyperball and it possesses quite a complicated structure. In the two 
previous sections we employed the presented procedure to a qubit (N = 2) and to qutrit 
(N = 3). The former is quite simple. On the other hand analysis of a qutrit gives 
support to all above given remarks. 

The employed representation of the density operator seems to be quite useful. We 
feel that we have shown that it can successfully be used to investigate issues which 
were not studied earlier and which seem to be of contemporary interest. Moreover, the 
procedure to investigate positivity is valid for any N. 

Finally, we would like to indicate some possibilities which probably deserve further 
attention. Polarization operators are spherical irreducible tensors jjj. Therefore, their 
tensor products preserve their character. This fact seems to be promising in the 
investigations of entangled states. We hope that this can provide new insights into the 
structure and geometry of entangled states. Moreover, as indicated by Biederharn and 
Louck jSj, polarization operators consitute just another (non-standard) set of SU(N) 
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generators. This set is endowed with an interesting stucture as it follows from careful 
inspection of the structure constants in (|20p. Investigation based upon this fact go 
beyond the scope of this work, but seem to be an interesting subject for future studies 
of the space of allowed generalized Bloch vectors and therefore of the geometry of the 
space of density operators. 

We end this paper with the hope that the revitalized expansion of the density 
operator in terms of polarization operators will prove useful in other investigations. The 
presented discussion of the correspondence between density operators and generalized 
Bloch vectors (especally in the light of requirement of positivity) is applicable to 
any dimension. The quantities appearing here are strongly connected with angular 
momentum theory and thereby well-known. This, in our minds, indicates that the 
representation discussed here may be more practical then the "standard" one. 
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